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We report results from a molecular dynamics (MD) simulation on the conformations of a long 
flexible polyelectrolyte complexed to a charged sphere, both negatively charged, in the presence of 
neutralizing counterions in the strong Coulomb coupling regime. The structure of this complex is 
very sensitive to the charge density of the polyelectrolyte. For a fully charged polyelectrolyte the 
polymer forms a dense two-dimensional "disk", whereas for a partially charged polyelectrolyte the 
monomers are spread over the colloidal surface. A mechanism involving the overcharging of the 
polyelectrolyte by counterions is proposed to explain the observed conformations. 
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Polyelectrolytes in polar solvents are polymers carry- 
ing dissociated ionic groups. When a polyelectrolyte is 
in the vicinity of a charged colloidal particle, both may 
coagulate leading to charge complexation. Studying this 
process is motivated by many sources. The presence of 
polyelectrolytes has important effects on the stabilization 
of colloidal suspensions Q|. Besides that, for polyelec- 
trolytes such as DNA the interaction with the interface 
of charged membranes or charged particles (histones) is 
crucial for many biophysical properties Finally long 
range Coulomb interactions represent a theoretical chal- 
lenge, especially for the understanding of effective attrac- 
tions between like charged bodies (|, || |[ |[. 

The adsorption of polyelectrolytes onto an oppositely 
charged spherical particle has recently been experimen- 
tally extensively studied 0, ||, 9|. Various authors 
have also investigated this phenomenon theoretically 
[10, O, [l2], [ll| |l4], [l5|, [l(| and by numerical simulations 
[17, [L8l |19|. However, much less is known concerning 
the complexation of a charged sphere with a like-charged 
polyelectrolyte. To our knowledge there has been no 
study in this direction until now. 

In this Letter, we report the rather unexpected com- 
plexation between a charged sphere and a long flexible 
polyelectrolyte, both like (here negatively) charged. This 
article constitutes a first attempt to elucidate this strik- 
ing phenomenon. We present results of MD simulation 
of the two macroions taking into account the counterions 
explicitly, but add for simplicity no salt. We propose a 
mechanism stemming from the polyelectrolyte overcharg- 
ing to explain the complexation structure as well as the 
observed polyelectrolyte conformations. 

The MD method employed here is based on the 
Langevin equation and is the same as the one employed 
in previous studies pC] . Consider within the framework 
of the primitive model one spherical macroion character- 
ized by a radius tq and a bare charge Q = —Zj^e (where 
e is the elementary charge and Zm > 0) surrounded by 
an implicit solvent of relative dielectric permittivity e r . 
The polymer chain is made up of N m monomers of diam- 
eter I. Both ends of the chain are always charged. The 



monomer charge fraction is / (i.e. every 1// monomer is 
charged) so that the chain contains N cm = (N m — 1)/+1 
charged monomers. The monomer charge is q m = —Z m e 
(with Z m > 0). To ensure global electroneutrality we 
added N c small counterions of diameter a and charge 
+Z c e (with Z c > 0). The whole system is confined in an 
impermeable spherical cell of radius R, and the spherical 
macroion is held fixed at the center of the cell. The dielec- 
tric permittivity is the same everywhere, in and outside 
the cell. 

Except for the spherical macroion all particles are mo- 
bile. Excluded volume interactions are introduced via a 
pure short-range repulsive Lennard-Jones potential given 
by 
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for r — ro < 2 1 / 6 er, and otherwise. We have set 
ro = 7a for the macroion-counterion interaction, and 
ro = otherwise. The length and energy simulation units 
are defined by a (counterion diameter) and clj, respec- 
tively. The closest center-center distance of the ions to 
the spherical macroion is a — ?*o + a = 8(7. The macroion 
volume fraction is defined as /m = (a/R) 3 and was fixed 
with R = 40ct to 8 x 10" 3 . 

The electrostatic interaction between any pair ij, 
where i and j denote either a macroion or a charged 
microion (counterion or charged monomer), reads 
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where 1b — e 2 /4ireoe r kBT is the Bjerrum length. To 
link this to experimental units and room temperature we 
denote €lj = ksT (T = 298 K). Being interested in the 
strong Coulomb coupling regime we choose the relative 
permittivity e r = 16, corresponding to Ib — 10er (with 
a = 3.57 A), divalent microions (Z m = Z c = 2) and 
Z M = 180. 
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FIG. 1: Typical equilibrium configurations of the colloid- 
polyelectrolyte complex for (a) run A (/ = 1), (b) run B 
(/ = 1/2) and (c) run C (/ = 1/3). Monomers are in white 
and counterions in red. 



The polyelectrolyte chain connectivity is modeled by 
using a standard anharmonic finitely extensible nonlinear 
clastic (FENE) potential in good solvent (see for example 

Ref. which reads U f en Eif) = —\hR$\tl 1 — jjs , 

where k is the spring constant chosen as lOOOfcsT/cr 2 
and i?o = 1-5(7. These values lead to an equilibrium 
bond length I = 0.8a. 

Four different parameter combinations, denoted by run 
A, B, C and D, were investigated and are summarized in 



TABLE I: Simulation parameters for the runs A, B, C and D. 
f denotes the charge fraction, N m the total number of chain 
monomers, N cm the number of charged chain monomers, and 
N c the total number of counter ions. 



Run 


A 


B 


C 


D 


1// 


1 


2 


3 


5 


N m 


256 


257 


256 


256 


N cm 


256 


129 


86 


52 


N c 


346 


219 


176 


142 



Table |[ Going from run A to D we decreased the charge 
fraction / from 1 to 0.2, and thereby decreased the linear 
charge density Xpe of the polyelectrolyte. The contour 
length of the chain is much larger than the colloidal parti- 
cle diameter (N m l / (2ro) w 14 times), so that in principle 
the chain can wrap around the colloidal particle several 
times. 

Figure |l| shows typical equilibrium configurations of 
the colloid-polyelectrolyte complex. Note that in all re- 
ported cases complexation occurs and the polyelectrolyte 
is completely adsorbed on the colloidal surface. How- 
ever, the structure of the resulting complexes depends 
strongly on /. For the fully charged polyelectrolyte case 
[see Fig. [|(a) with / = 1] the chain monomers are closely 
packed together with their counterions forming a two- 
dimensional dense aggregate. This conformation consists 
of alternating lines made of monomers and counterions, 
respectively. When the linear charge density is reduced 
[see Fig. |](b) an d Fig. 0(c)], the complex structures are 
qualitatively different. In these cases the chain monomers 
are no longer densely packed. For run B [Fig. 0(b)], the 
monomers spread more over the colloidal surface and the 
polymer chain partially wraps around the sphere suggest- 
ing almost a surface pearl necklace structure. For run C 
[Fig. 0(c)], the chain wraps entirely over the colloidal 
surface, leading to a quasi isotropic distribution of the 
monomers around the spherical macroion. 

The physical reason of the complex formation is due 
to the strong counterion mediated attractions. Basi- 
cally, the charged species is trying to crystallize in a 
way which is compatible with the topological constraints 
(here mainly the chain connectivity and the macroion 
surface). 

To quantify the adsorption of the chain monomers and 
counterions on the macroion surface, we look at the ob- 
servable P(r), being defined as the amount of monomers 
(counterions) reduced to the total number of monomers 
N m (the total number of counterions N c ) within a dis- 
tance r from the spherical macroion center. Results are 
depicted in Fig. ^, where we observe for all runs that 
the particles are condensed within a distance of about 
lOtr from the colloid center and more than 80% of the 
monomers and counterions are within a distance of 9.3a 
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FIG. 2: Fraction of counterions (thin lines) and monomers 
(thick lines) adsorbed on the spherical macroion as a function 
of the distance from its center. 



FIG. 3: Radius of gyration of the polymer as a function of of 
the monomer charge fraction / (runs A-D). 



Xpe is defined as 



from the colloid center, corresponding roughly to two 
particle layers. Due to strong electrostatic attraction be- 
tween the sphere and the counterions and strong elec- 
trostatic repulsion between the sphere and the charged 
monomers, the first layer (r ~ a — 8<r) is exclusively 
made up of counterions. Note that the monomer deple- 
tion in this first layer also concerns neutral monomers 
(runs B-D) and this effect is attributed to the chain 
connectivity. The charged monomers in the second layer 
repel each other strongly, and keep for high values of / 
the chain under tension and thus the neutral monomers 
almost out of the first layer. The second layer consists of 
both, chain monomers and counterions. 

Next, we investigate the radius of gyration R g (in three 
dimensions) of the chain in order to gain insight of the 
spreading of the monomers over the sphere. Results are 
reported in Fig. || It clearly indicates that R g increases 
with increasing /, which demonstrates that the spreading 
of the monomers over the macroion surface is enhanced 
by decreasing the polyelectrolyte charge density. The 
jump in R g is particularly large between the results for 
/ = 1 and f < I. This is in agreement with the visual 
inspections of the polymer conformations presented in 
Fig. [l| Moreover, the isotropic case (monomers fully 
spread over the particle) corresponding to R g ~ a + I = 
8.8(7 is already reached for / = 1/2 (run B). 

In the following we are going to argue that the poly- 
electrolyte overcharging is a fundamental key to explain 
the observed complex structure. Let N c d be the number 
of condensed counterions on the polyelectrolyte where 
counterions are assumed to be condensed onto the poly- 
electrolyte when they lie within a distance r c = 1.2a to 
the nearest chain monomer. Then the overcharging ratio 
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which is merely the ratio between the amount of the to- 
tal condensed counterion charge (on the chain) and the 
polyelectrolyte bare charge. 

This overcharging can also be analytically predicted by 
the simple assumption that the uncondensed counterions, 
i.e. those who are not "attached" to the chain, have the 
same surface pair distribution as in the case where the 
polyelectrolyte and its condensed counterions are absent, 
i. e., where we just have the macroion and its condensed 
counterions present p3| . This allows us to consider the 
polyelectrolyte chain as a neutral object, namely the bare 
charged chain plus its own neutralizing counterions which 
then in turn gets overcharged by intercepting all counteri- 
ons of the macroion which are thought of being uniformly 
distributed over the macroion surface. We consider coun- 
terions as belonging to the chain if their center lies within 
a ribbon of width 2r c and area A r n, = 2r c N m l. If c is 
the counterion (of the macroion) concentration then the 
theoretical overcharge xth is merely given by 
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and for N m > 1 we have 



Xth ~ 1 + C/f, 



(4) 



(5) 



where C 

Results are presented in Fig. H. It indicates that in 
all cases overcharging is present (i.e., x.pe > 1), and 
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FIG. 4: Polyelectrolyte overcharge as a function of the 
monomer charge fraction / (runs A-D). The dashed line cor- 
responds to the theoretical prediction where Eq. wm was used. 



that xpe increases with decreasing polyelectrolyte charge 
density. We have excellent agreement (less than 10% 
difference) between the simulational results and our toy 
model [Eq. (j2|)]. In turn it explains why xpe varies quasi 
linearly with 1//. 

The /-dependency of the complexation structure can 
be explained with this overcharging. Indeed, in terms 
of many-body physics the overcharged polyelectrolyte can 
be seen as a dressed (or renormalized) particle [bare par- 
ticle + clothing] with a new effective (or renormalized) 
linear charge density X* PE = —(xpe — l)Ap_E which has 
an opposite sign to the bare linear charge density Xpe 
[p3. In this respect, one can also define a renormalized 
charged monomer of effective charge 
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Using Eq. (||) and the results of Fig. [| this shows 
that the absolute value of increases with increasing 
1//, and therefore the renormalized polyelectrolyte self- 
interaction leads to a stronger effective monomer- effec- 
tive monomer repulsion, which in turn explains why the 
chain expands with increasing 1// (see Fig. gand Fig. | 
for the corresponding structures) p5| . 

In summary, we have shown that in the strong 
Coulomb coupling regime, a charged spherical macroion 
and a long polyelectrolyte can make a complex even if 
both carry many like charges. The resulting structure 
of the charge complex depends strongly on the linear 
charge density of the chain. For a fully charged poly- 
electrolyte, its chain monomers and condensed countcri- 
ons are densely packed and the polymer wraps only par- 
tially the sphere. By decreasing the linear charge density, 
the wrapping becomes complete. The resulting polyelec- 
trolyte conformations can be explained by the degree of 
overcharging of the polyelectrolyte which in turn depends 
on its linear charge density. 



A future study will include many other important ef- 
fects, such as chain length and flexibility, enlarged per- 
mittivity constant (weaker Coulomb coupling), added 
salts, and microions valence to name just the most impor- 
tant parameters. Nevertheless, our observations should 
trigger new theoretical and experimental works. 
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